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Abstract 

QCD corrections to the QED formula for parton distribution functions of the longitu- 
dinal virtual photon are derived in the leading-logarithmic approximation. It is shown 
. that the resulting PDF satisfy the same homogeneous evolution equations as those of 

\ hadrons, but contrary to the latter are perturbatively calculable. Properties of these 

distribution functions are discussed and their phenomenological consequences outlined. 
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O ' 1 Introduction 

<N . 

Parton distribution functions (PDF) of the photon can be written as sums of the hadronic (HAD) 
' and pointlike [] (PL) parts. The former satisfy the same homogeneous evolution equations as those of 

hadrons, whereas the latter satisfy the inhomogeneous ones. Within a subset of pointlike solutions 
! specified by the value of the scale Mq at which they vanish we can thus write 

D(x,M ,M) = D HAD {x,M ,M) + D PL {x,M ,M), D = q,q,G. (1) 

O ■ 

Due to the freedom in the choice of Mo the separation of PDF into their pointlike and hadronic parts 
is, however, ambiguous. In [jl] we discussed practical aspects of this ambiguity for the Schuler- 
Sjostrand sets of parameterizations 0. The quark distribution functions of the photon can be 
Q_i' written as combinations of singlet and nonsinglet parts 

<D : 
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Due to the presence of hadronic contributions, individual qf s are in general independent nonsinglet 
distributions. Their pointlike parts are, however, proportional to the corresponding charges and 
(neglecting quark mass effects) pointlike parts of qi can therefore be expressed in terms of one 
singlet and one nonsinglet quantity 

l! L = i£ - (e 2 )) Q NS ' PL + <eV' PL , S = 2n / <eV< PL (3) 

The pointlike part of g NS is defined via the resummation of multiple gluon emissions off the quark 
from the primary 7 — ► qq splitting. In units of a/27r 
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1 Frequently also called anomalous, because they contain the term that is absent in the case of hadrons. I prefer 
the denomination pointlike, which reflects the fact that they arise from pointlike 7 — > qq coupling. 
2 Throughout this paper the expressions for PDF of the photon will be given in these units. 
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where /cns = 3(x 2 + (1 



2 ) and a s (M) satisfies the standard renormalization group equation 



IKF =/3(« S (M)) 
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where, in QCD with nj massless quark flavours, the first two coefficients, (3$ = 11 — 2nj/3 and 
/?i = 102 — 38rij/3, are unique, while all the higher order ones are ambiguous. In this paper only 
the LO QCD corrections will be considered and therefore only the first term in (||) will be taken 
into account. In terms of standard moments the series in (0) resums to 



,NS,PL, 



n,M ,M) 



a s (M) 
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a s (M )J 
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For transverse virtual photon jt, the above considerations can be repeated with only a minor 
modification reflecting the fact that its momentum P 7 is off-shell, P 2 = —P 2 > 0. The exact 
meaning of the quantity r denoting in (||) the outgoing parton off-shellness is somewhat ambiguous. 



stand for parton momentum and mass 



I have taken r as r = — (p 2 — m 2 )/x, where p and m L 
respectively. Working out the kinematics of the 7* — > qq vertex in the collinear limit yields r r 
P 2 + m 2 /(x(l — x)). For massless quarks this implies replacing in 



Ml with P 2 . 



2 The concept of PDF of 7^ in QED 



Before discussing PDF of 7£ in the framework of QCD, let me illustrate the usefulness of this 



concept on the QED analysis of DIS of electrons on virtual photons. In the region m 2 <C P 2 <C Q 2 



experiments measure the following quantity 
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where the cross-sections ajk are in general functions of W 2 , P 2 and Q 2 and : 
The 0(a) QED expressions for individual contributions to read P] 



Q 2 /(W 2 +Q 2 +P 2 



F^j, (x, P ,Q ) 

F^ D (x,P 2 ,Q 2 ) 
F^ D (x,P 2 ,Q 2 ) 



Q 
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V v ' 
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F L Q L ED (*,P^ 



nonpartonic 



(8) 

0) 
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The denomination "partonic" refers to the fact that the corresponding term arises from integration 



over r of the quark from the primary splitting 7 — ► qq close to its lower limit T n 



P\ i.e. 



from the region of phase space where the emitted quarks (antiquarks) are almost collinear with 
the incoming photon. In the case of F^P part of the integrand is proportional to 1/r, which 
yields the logarithmic term, whereas the "partonic" 2x(l — x) — 1 = —(x 2 + (1 — x) 2 ) in F^ D and 
4x(l — x) in F^l^ arise from integration over P 2 /t 2 . In both cases the resulting expressions have 



The first and second indices corresponding to probe and target photon respectively. 
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a clear parton model interpretation. On the other hand, the "nonpartonic" contributions result 
from integration over the whole phase space and cannot be interpreted in terms of PDF. Adding 
all terms in ([?]) we get 



F<f D (*,p 2 ,Q 2 ) = 6 



Q 



2 



{x 2 + (l-x) 2 ) ln^ + 12x(l 



(11) 



One might argue that there is no reason to distinguish the origins and interpretation of the indi- 



vidual terms that add up to the nonlogarithmic part of (11), and consequently, the effects of target 
72 can be included as part of the "constant" term in (11). Although legitimate, I prefer not to 
adopt this procedure and keep track of the origins and interpretation of individual contributions 
to the constant part of Fj s because: 

• The "partonic" terms, coming from the vicinity of mass singularities are universal quantities, 
characterizing the incoming 7^, whereas the "nonpartonic" ones are process dependent. 

• Virtuality dependence of different terms is different. The above expressions hold for m 2 <C 
P 2 <C Q 2 - As P 2 /m 2 — > the nonpartonic terms, coming from integration over the whole 
phase space, remain unchanged, whereas F^ L vanishes by gauge invariance as Fj, L oc P 2 /m 2 
In the case of D the partonic part of the constant term gets additional contribution from 
the term m 2 /(m 2 + x(l — x)P 2 ) — * 1, absent for P 2 ^> m 2 . Taking these facts into account 

allows us to derive in a physically transparent way the constant term of i^ ED (x, Q 2 ) for the 
real photon: [2x(l - x) - 1] + [2a; (1 - x) - 1] + 4x(l -x) + l = 8x(l - x) - 1. 

The fact that F^° = 24x(l — x), which does have a parton model interpretation coincides with 
F® T , which does not, reflects the symmetry of the total cross section gtl{i* {Q 2 )l* (P 2 )) with 
respect to P 2 and Q 2 . 



3 Pointlike part of 7^ in QCD — nonsinglet distribution 

The purely QED prediction for the pointlike part of q^ s depends on the quark mass which 
provides the scale governing its threshold behavior. In this paper I consider QCD corrections to 
QED formula for massless quarks (or for x(l — x)P 2 ^> m 2 ) and in the region P 2 M 2 . For 7£ 
the resummation of diagrams in Fig. |l| in the collinear region takes the form 

p 2 ,m 2 ) m M „ /; anf + £ 4* (=) i: £ anf 

where ki(x) = 12x(l — x) and the first term in (^) defines the QED contribution 

«r<^.M 2 ) - "lW (~^)= (l - 5) P -C (13) 

Performing in (12) the inner integrals over t\ we get 

4 In not stated otherwise all expressions for quark and gluon distribution functions mentioned in the following 
concern their pointlike part and the superscript "PL" will therefore be dropped. 
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Figure 1: Diagrams denning the pointlike part of nonsinglet quark distribution function of the 
photon. The resummation involves integration over quark virtualities r < M 2 . 



qf s (x,P 2 ,M) = k L (x)(l 



p2 

M 2 
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1 dw (0) /x 



1 dy p (0) 
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p2 r 2 2ir 
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r 2 27r 
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Note that (1-P 2 /M 2 ) and (l-P 2 /r 2 ) are defined also outside the region P 2 < M 2 (or P 2 < t 2 ), but 
are negative there. In order to avoid this unphysical region, all integrals in ( |l4| ) will be understood 
to be multiplied by 9(1 — P 2 /M 2 ) (or 9(1 — P 2 /r 2 )), guaranteeing thus the vanishing of ( |T4l ) as 
P 2 — > M 2 . Taking into account in each brackets of (|TJ) only the unity, we get 



q r( X ,p 2 ,M 2 ) 



k L (x) 



poJx y \yj 



1(1 



1 ^p(o) ^ 
u; qq \w 



1 ^p(O) 
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where s = ln(M 2 /A 2 )/ln(P 2 /A 2 ). In momentum space (15) resums to 



^ s (n,P 2 ,M 2 ) 



k L (n) 



a s (M 2 



2P, ( ° ) (n)//3 



a s (P 2 ) 

which satisfies standard hadronic evolution equation in the nonsinglet channel 



dgg a (n,P 2 ,M 2 
dlnM 2 



^^P£Kn) q T(n,P 2 ,M 2 



(16) 



(17) 



The resummation in (14) resembles the definition of nonsinglet quark distribution function of 
hadrons, with Ul(x) playing the role of "bare" quark distribution function and P 2 the initial scale. 



Note, however, that (|iq ) has been derived under the assumption P <C M . Translating (|iq ) into 
x-space by means of inverse Mellin transformation leads to scaling violations illustrated by solid 
curves in Fig. ^. Multiple gluon emission leads to softening of the distribution with increasing M 2 , 
a feature typical for hadrons. 

As emphasized above, the result (|16D holds for P 2 <^ M 2 only. For P 2 — ► M 2 the expression 



(14) vanishes by definition, whereas (|l^) goes to kj_,. Anticipating the results of the next paragraph 
let me approximate the inclusions of the terms neglected in ( |l6| ) by replacing the first term on the 
r.h.s. of (|l5|) by ki(x)(l — P 2 /M 2 ). This amounts to replacing (|l6|) by 



NS ^n,P 2 ,M 2 ) 



k L (n) 



a s (M 2 ) 
a s (P 2 ) 



2P£ ) (n)/(3 



(18) 



and leads to series of dotted curves in Fig. |2[ As expected these curves are close to the solid ones 
for P 2 <^ M 2 , but substantially different for P 2 approaching M 2 from below. As emphasized at 
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Figure 2: Left: Scaling violations for q^ s . Solid curves correspond to ( |T6| ) for P 2 = 2 GeV 2 and 
(from above at large x) M 2 = 3, 10, 10 2 , 10 3 , 10 4 GeV 2 , the dotted ones to the same values of M 
but the formula (18). For the highest two values of M 2 the dotted curves are indistinguishable from 
the solid ones. Right: Comparison of the pointlike part of as given by ( |l6l ) with the hadronic 



part of q^ s evaluated from SAS2D parameterization of ut according to the formula (| 



the beginning, PDF of virtual photons make good sense for P 2 <C M 2 only and Fig. ^ suggests 
that this means, roughly, P 2 < 0.2M 2 . On the other hand, from the point of view of mathematics 
the evolution equation (17) as well as its solution fli~6| ) can be considered for arbitrary P 2 and M 2 . 
In analyzing higher order corrections to ( |l6|) I will make use of the following approximation 



V(x ) 



dx 1 



1 



1 



, Xo x 2 ln(x) x ln(2x )' 
which for 2 < xq <C xm is very accurate and allows us to write 
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/3o J P 2 r 2 ln(r/A 2 ) 
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x 



2 / a2 



x 2 ln(x) Po ln(2xo) 



2tt 
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+ 0(a 2 (P 2 )) (20) 



P 2 /A 2 , x M = M 2 /A 2 . This integral first appears in evaluating the 



contribution of the second part of the 0(a s ) term in (|l4|), which has been left out in deriving (^ 

2 dr 2 a s (T 2 ) . a s (P 2 



1 ^p(O) ( f 



F 2 



In momentum space (21) turns into 
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2tt 
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a s (M 2 ) 
a s (P 2 ) 



-,-2P 9 ( 9 0) (n)//3 



(21) 
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where the expression after the arrow follows when in all terms on the r.h.s. of ([l4|) the inner integral 
over dr is approximated using ( ^9|) as (a s (P 2 ) / P 2 — ol s {ti)/ti) and only the first term retained. 
Starting from the third term on the r.h.s. of ( |l4l) the second term, i.e. cx s (ti)/ti is not negligible as 
the upper limit t\ does not coincide with M 2 , but spans the whole range P 2 < ti < M 2 . Taking it 
also into account yields (]i~6| ) multiplied by a term proportional to a 2 (P 2 ). Repeating this procedure 
at higher orders leads to the general expression 



q™(n,P',M') = k L {n) \ 1 + >>(n)aJ(P 



k=l 



a s (M 2 ) 
a s (P 2 ) 



2P™(n)//3 



(23) 
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G PL (x,M 2 ) 

Figure 3: Part of the set of diagrams defining the gluon distribution function of the photon. 

where c\ = —P^[n)/2ir etc. Although calculable all the coefficients Ck(n),k > 1 were mentioned 
only to indicate the origin of the general structure of the results after the integrations in (|l4|) , but 
should be discarded in the leading-log approximation, used in deriving (|l6|). 

To estimate the hadronic part of q^ s , I have assumed that its onset is governed by the same 
ratio P 2 /niy of the photon virtuality P 2 and vector meson mass rriy as the decrease of q^ s with 
increasing P 2 : 

qf < HAD (z,P 2 ,M 2 ) ^ r(PV S ' HAD (*,P 2 ,M 2 ), r(P 2 ) ^ (pgq^r) ■ (24) 

In Fig. |2|b (|l6|) is compared to (|24|) assuming <^ S ' HAD = xu(x)/3e 2 = 3xu(x) and using the SAS1D 
parameterization of u(x,P 2 , M 2 ) pi Estimated in such a way, the hadronic part of seems to be 
safely negligible for P 2 > 1 — 2 GeV 2 . 



4 Pointlike part of 7^ in QCD — singlet and gluon distributions 



The fact that satisfies homogeneous evolution equation suggests the same for the quark singlet 
and gluon distribution functions. Indeed, following the procedure of previous Section one can resum 
the series coming from diagrams in Fig. [3| as follows 



G L (n,P 2 ,M 2 



2n f (e 2 )k L (n)P^(n 



a s (M 2 
a s (P 2 ) 



2M -2PW(n)/A> 



2n f (e 2 )G s 1 



(25) 



Adding multiple gluon emissions on the outgoing quark line modifies (25) as follows 

"a s (M 2 )l- 2P ^ (n)//3 ° \a s (M 2 



Gl(n,P',M< 



k L (n)P§}(n) 



P 



Gq_ 
P 



GG (n)-P^{n) 

Combined with q^ s defined in (|l^) G|(n,P'M 2 ) satisfies the evolution equation 
dG s L (n,P 2 ,M 2 ) _ a s {M 



a s (P 2 ) 



a s (P 2 ) 



(26) 



dlnM 2 



2tt 



Pi°l(n)Gl(n, P\ M 2 ) + P$ P\ M 



(27) 



As for q^ s ( |26|) was derived assuming P 2 « M 2 , but from the point of view of mathematics 
can be considered for arbitrary values of P 2 < M 2 . Note that whereas q^ s (n, P 2 , P 2 ) = ki(n), 
G^(n, P 2 , M 2 ) defined in (26) vanishes for P 2 = M 2 . The equation (^?j) is not exactly the evolution 
equation for gluons in hadrons as the latter contains quark singlet instead, but the diagrams 

5 In this way the hadronic part of qr S is in fact overestimated. For P 2 > 2 GeV 2 the suppression factor P 2 /{P 2 + 
my) is close to unity and has therefore only small effect. 
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Figure 4: Scale dependence of g| and G| for nj = 4 and A 2 = 0.1 GeV 2 . Ordered from below at 
small x the curves correspond to M 2 = 3, 10, 10 2 , 10 3 and 10 4 GeV 2 respectively. 



included in the definition ( f26|) do not exhaust all possible ladder diagrams defining The above 
considerations are, however, sufficient to conclude that "dressed" PDF of are constructed from 
the "bare" ones in the same way as those of hadrons and therefore satisfy the same evolution 
equations as the latter, with initial conditions at M 2 = P 2 having the form 



qf S (x,P 2 ,P 2 ) = ql(x) = k L (x) = 12x(l - x), G S L (x,P 2 ,P 2 ) 
Using the formulae of M we can write explicitly 
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Kin) 



a[n) 



d±(m) 



yj (Pff(n) - 4gH) 2 + 8n f P^(n)P^(n), 



P$\n) - P®{n) + K{n) 
2K(n) 



> e ( n ) 



qG 

Gq' 



?$q\n) 



-P 3 ( ?(n, 



P^(n) ± K (n) 



1% 



(31) 
(32) 
(33) 



Together with qf^ as given in ( [l6|) this completes the evaluation of PDF of "f* L . The singlet charge 
factor 2n/(e 2 ) factorizes out of both and Gl due to vanishing of Gl at the initial scale P 2 . The 
nontrivial dependence on nj remains inside n{n). 

Performing numerically the inverse Mellin transformation of ( ^9| - |30| ) leads the scale dependence 
of gf and G L in x-space, illustrated in Fig. || for P 2 = 2 GeV 2 , A 2 = 0.1 GeV 2 and rif = 4. 
The resulting parameterization [] of PDF of j L based on eqs. ((H), (H) and (H) contains A as a 
free parameter. The fact that PDF of satisfy the homogeneous evolution equations implies the 
validity of standard momentum sum rule. Recall that its violation for ^ is due to the fact that 
the "bare" quark distribution function of 7^, g^ ED oc InM 2 , depends on the factorization scale M. 
6 Obtainable upon request from chyla@fzu.cz. 
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Figure 5: Comparison of the contributions Ftt and Ftl of target 7^ and to F e g (left) and D e g 
(right) for P 2 = 2, 5 GeV 2 and A/ 2 = 100 GeV 2 as given by QED and QCD formulae flU) and @ 
for 7£ and SAS1D parameterization for 7^. 

The QCD induced scale dependence of PDF of 7^ is formally reminiscent of the idea of "dy- 
namically generated" PDF of hadrons advocated by the GRV group (see || and references therein). 
Note, however, that the initial conditions ((2^) themselves have no direct physical interpretation as 
all our formulae were derived under the assumption P 2 <C M 2 . 



5 Phenomenological consequences 

In H we discussed the numerical relevance of the contribution of resolved 7^ for two physical 
quantities: F^(x,P 2 ,M 2 ) and the effective PDF 

n 1 

D eS (x, P 2 ,M 2 ) = £ P 2 ,M 2 ) + q t (x, P 2 ,M 2 )) + |g(x, P\ M 2 ) (34) 
i=l 

relevant for dijet production in ep and e + e~ collisions. The contributions of 7£ plotted in Figs. 
2 and 3 of and reproduced by dashed curves in Figs. [5] below were calculated using simple 
QED formula (13) for the quark content of 7?. They peaked around x ~ 0.7 and were more 



pronounced for F^ than for -D e ff- QCD corrections to PDF of 7^, given in eqs. (16),(p9|) and ( 
modify this situation significantly. Compared to purely QED predictions, the QCD effects suppress 
quark distribution functions of 7^ at large x and enhance them on the other hand for x < 0.4. 
This enhancement with respect to the contribution of 7^ is in fact so large that for x < 0.02 the 
contribution of to i 7 ^ exceeds that of 7^. For D e d the region x < 0.4 is dominated by gluons 
from both 7^ and 7^, and again the contributions of 72 is comparable to those of 7^. The relevance 
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of 72 is further enhanced with increasing P 2 (but keeping P 2 <C M 2 ). but decrease with increasing 
M 2 . Taking into account experimental conditions at LEP and HERA, Figure [5] suggests that dijet 
production at HERA ]jj offers particularly suitable place to look for the effects of 7^ . The quantity 
D e g allows us to express the relevant cross sections in a simple but only approximate way. Work 
on incorporating QCD improved PDF of 7^ within the LO Monte-Carlo event generators as well 
as NLO parton level calculations is in progress. 
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